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A full, high-statistic 7 = and 7 = 2 n — n scattering calculation including the disconnected 
diagram is performed on 2+1 flavor, domain wall fermion gauge configurations generated using 
the Iwasaki gauge action with j3 = 2.13, a 16 3 x 32 volume with L s = 16 and aT x = \.13GeV . 
Using the same light quark propagators and additional strange quark propagators, we study the 
77 and tj' mesons, where disconnected diagrams also make an important contribution. Our % — 
K calculation shows a good exponentially decaying signal from the disconnected graph, which 
makes a small contribution to the 1 = scattering amplitude at least for the short times where it 
can be accurately computed. We are able to resolve the rj and rj' states and see the pattern of 
SU(3) flavor symmetry breaking found in Nature. 
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1. Introduction 

Disconnected diagrams play a very important role in many lattice QCD calculations, such as 
the calculation of 7], 77 and a masses, the isospin % — % scattering length, e' from K to Tin 
AI = 1/2 process, etc.. Since the correlation functions corresponding to these diagrams reflect the 
quantum fluctuations of QCD vacuum, the noise does not exponentially decrease with increasing 
time separation between source and sink. Thus, the signal present in these correlation functions 
(which behaves like exp (—m e fft) for large time separations) is quickly buried in the noise and it is 
very difficult to get useful information from large time separation. Therefore, a good understanding 
of the disconnected diagrams requires a full dynamical calculation and large statistics. Several 
attempts related to these calculation have been done for the r]' in the case of SU(2) [1, 2, 3] and 
the a [4, 5]. It would be more interesting to calculate the physical T] and r\' masses in the flavor 
symmetry broken SU(3) case where m s > m/. 

The n — 71 scattering length calculation per se, is very important to expand our understanding 
of the strong interaction on lattice QCD. It also represents a first step toward a more interesting 
calculation of the direct CP violation measure e' which is believed to be a quantity that may call for 
new physics. While much work has been dedicated to the 1 = 2 channel of % — % scattering [6, 7], 
few lattice QCD calculations have been done for the 1 = case, possibly because the disconnected 
diagram makes the problem much harder than 1 = 2. Y. Kuramashi et. al. studied the isospin 
channel % — % scattering [8] almost two decades ago, but he ignored the disconnect diagram and 
also used the quenched approximation. In this paper, we calculate all the diagrams, especially 
providing a close look at the disconnect diagram. We will get a rough idea about how large the 
disconnected diagram's contribution is. 

2. Details of Lattice Calculation 

Our calculation is a full unitary calculation and is based on the RBC/UKQCD 16 3 x 32, 
L s = 16, 2+1 flavor domain wall fermion, j8 = 2.13 Iwasaki gauge action lattices. The inverse 
lattice spacing for these lattices is determined to 1.73(3)GeV [9]. We have three ensembles of 
such configurations with light sea quark mass wi/=0.01, 0.02, 0.03 respectively. They all have the 
same strange sea quark mass m s = 0.04. The ensembles with mi = 0.01 and 0.02 each has 150 
configurations, and the m/ = 0.03 ensemble has 281 configurations. 

We use a wall source and a wall sink for the propagator. The calculation of the correlation 
function for the disconnected diagrams naturally requires us to calculate the propagators on the 
time slices of both source and sink, so we decide to calculate the propagators on all time slices, 
which requires the calculation of 32 separate propagators in our case. By doing so, the statistics 
are greatly improved since we can put the source at all possible time slices. In a sense, we extract 
all possible information from a given configuration. Finally, we exploit the need to compute many 
propagators on each configuration by using Ran Zhou's eigenvector accelerator code, speeding up 
the computation of the propagators by 60 percent. 

3. The r\ and r\t masses 

Using approximate SU(3) flavor symmetry, we expect that we can create the rj with the octet 
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operator 0g, and the r\' with the singlet operator 0\, 

j] «0 8 = -=(uy$u + dy$d-2sy$s) (3.1) 
v6 

T}' k.O\= —(uysu + dYsd + sYss) (3.2) 
v3 

We will try to demonstrate this is true from lattice QCD calculation. For future convenience, we 
define a light operator and a strange operator, 

Oi = -^=(iiY5ii + dYsd) (3.3) 

O s = sy 5 s (3.4) 

and consider the correlation function matrix of these two operators, 

( (0/(00/(0)) (0,(00/(0)) \ = / Q-2D, -V2D X \ 
{ > { (0/(00.1(0)) (O s (0Oj(0)) J { -V2D X C s -D s J 

where Q, C s , D/, D s , and D x stand for all the five possible contractions which are shown in Figure 1. 
Then by diagonalizing the matrix C~ l (to)C(t) with a reference time to = 1, we expect to get the 
correlation function for 77 & 77 ' states, while the rotation matrix that diagonalizes this matrix gives 
the mixing angle for the light and the strange parts, 

77 = cos dri-Oi- sin 6^ ■ O s (3.6) 
7]' = sin 6 n > ■ Oi + cos d v > ■ O s (3.7) 

In the SU(3) flavor symmetry limit, sinfyj ^ = y| = 0.8165, cosfy ^ = = 0.5774. 

The effective mass plateau of 77 and v\' for m/ = 0.01 ensemble are shown in Figure 2. Using 
the range 3 < t < 7 where there are good plateaus, we get = 0.397(17) = 687(29)MeV, and 
m^i = 0.654(80) = 1.13(14)GeV. Figure 3(a) shows the mixing angle calculated by diagonalizing 
the correlation matrix. The data points with small time separation where we have good signal for 
the rj' agree with the SU(3) symmetry limit. This demonstrates that the T] is created by the octet 
operator and the 77 ' by the singlet operator, approximately. This ability of a full QCD calculation 
to reproduce the anomaly dominated SU(3) symmetry breaking pattern of the physical r\ and t]' is 
an important test of lattice QCD. 

Since the SU(3) flavor symmetry is approximately true, using Eqn. 3.1 and Eqn. 3.2, we can 
calculate rj and T]' correlation function directly from the 0$ and 0\ operator, 

(T 7 (/)T 7 t (0)) = i(Q + 2C,)-|(D,+D,-2D,) (3.8) 

(T 7 / (/)T J /t (0)) = ^(2Q+C s )- l -(4D I + D s + 4D x ) (3.9) 

this gives m n = 0.392(12) = 678(20)M<?V, = 0.671(82) = 1.16(14)GeV, and they agree with 
the results of correlation matrix diagonalization method. 

Our result for the 77 mass is consistent with Gell-Mann-Okubo formula, which states that 
3m^ +ml = 4m 2 k , and gives m n = 0.3855(25) = 667(4)MeV , where we have used m n = 0.2472(10) 
All.lMeV and m k = 0.3560(20) = 615.9M<?V. 
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Figure 1: The left panel shows all five diagrams that contribute to the r\ and r/' correlation functions. They 
are Q, C s , D/, D s , and D r respectively from the top left to bottom right. The solid line with an arrow stands 
for a propagator, and the black solid circle stands for a 75 insertion. The right panel shows our result for 
these diagrams for the m/ = 0.01 ensemble. 
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Figure 2: Left panel: effective mass plateau for the rj and the mass obtained using the fitting range [5:14]. 
Right panel: effective mass plateau for the r/' and the mass obtained using the fitting range [3:7]. 



For the v\' mass, since we can only get a reasonable signal from short time separation where 
the excited states may exist, our result may not be fully reliable. Figure 3(b) shows m„> calculated 
for mi = 0.01 and 0.02 ensembles. By assuming a linear relationship between and m n 2, it 
seems our result does agree with the experiment physical tj' mass. 



4. Two-Pion Scattering 

The diagrams that need to be calculated in order to study % — % scattering in both 1 = 2 and 
1 = channels are shown in the left panel of Figure 4, identified as the Direct(D), Cross(C), Rect- 
angular^), and Vacuum(V) diagrams. The right panel shows the actual data for these four con- 
tractions for the mi = 0.01 ensemble. We can see that diagram D makes the biggest contribution, 
then diagram C and R, and diagram V makes the smallest contribution. It is extremely noisy for the 
disconnected diagram(V), but still we can get a clear signal up to time separation 6. By increasing 
the statistics, or reducing the pion mass(making the 1 = signal fall less rapidly), we could obtain 
a better result for the V diagram. 
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Figure 3: (a). The Mixing angle, which determines the two eigenvectors of the matrix C _1 (fo)C(f). The two 
horizontal lines give the value of cos(0) and sin(0) in the SU(3) flavor symmetry limit, (b). The r\' mass 
from the mi = 0.01 and mi = 0.02 ensembles. The physical mass 957MeV is also labeled on the figure. 
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Figure 4: The left panel shows the Direct(D), Cross(C), Rectangular(R), and Vacuum(V) diagrams from the 
left top to right bottom. The right panel shows our result for these diagrams for the m; = 0.01 ensemble. 



The K — 71 correlation function for / = 2 and 1 = can be expressed in terms of these diagrams, 

< h{t) |/ 2 (0) > = 2(D-C) (4.1) 
</ (f)|/o(0) > = 2D + C-6R + 3V (4.2) 

We can get the energy of / = 2 and / = % — 71 state from the plateau of the effective mass plot, and 
the results are shown in Table 1. Since the disconnected diagram makes only a small contribution 
to the correlation function of / = scattering state but contributes very big noise, leaving it out 
makes the numerical result more precise but in trade with an unknown systematic error, and it is 
listed as (/ = Vout) in the table. 

To get the scattering length, we use Liischer formula [10, 1 1], which relates the energy of the 
two particle states with the scattering length, 

A£ = £-2^ = - 47r ^[l+c 1 ^ +c 2 (^ , ) 2 )]+O(l/L 6 ) (4.3) 

T)%-jrLi Li Li 



5 



Disconnected diagrams 



Qi Liu 



Table 1: The energy difference between the %—% isospin eigenstates and 2m K . The last column shows 
Ae(I — 0) when we leave the disconnected diagram out. The numbers in the bracket [] give the fitting range. 



ni[(conf) 


m n /[5 : 14] 


A £ (/ = 2)/[5:14] 


A £ (/ = 0)/[5:8] 


A £ (/ = Vout)/[5:12] 


0.01(150) 


0.2472(10) 


0.0197(14) 


-0.041(80) 


-0.0439(72) 


0.02(150) 


0.3248(9) 


0.0158(7) 


-0.086(79) 


-0.0268(69) 


0.03(281) 


0.3895(6) 


0.0135(5) 


0.04(10) 


-0.0155(50) 



Table 2: The K — K scattering lengths for the 7 = 2 and 7 = channels. The data in the bracket [] shows the 
minimum and maximum value corresponding to one standard deviation. The last column shows the result 
for I = scattering length when the disconnected diagram is omitted. 





ao{I = 2) 


a (I = 0) 


a Q (I = Vout) 


0.2472(10) 


-1.26(7) 


4.6 [-2.1:8.1] 


4.9 [4.2:5.4] 


0.3248(9) 


-1.31(5) 


7.9 [0.86:10.1] 


4.1 [3.0:4.9] 


0.3895(6) 


-1.34(5) 


-2.9 [-6.0:7.3] 


2.8 [1.7:3.8] 



where ci=-2.8373, C2=6.3752. The calculated scattering lengths are shown in Table 2 and Figure 5. 

Let us compare our result with theory preditions. The tree level #PT formula for the 71 — % 
scattering length is 

1 1nv\ , , 1 7m~ 

which we evaluate using the results for f % and m n obtained in Ref [12], 

f % = 0.0765 + 1 .02(/n z + 0.00308) , m\ = 2 x 2.285(m, + 0.00308) (4.5) 

The #PT results are shown in Figure 5 as the solid line, in comparison with our lattice QCD 
results. Notice that the tree level formula agrees well with the experiment and our results for the 
1 = 2 scattering length, while giving a value for I = around 30 percent smaller compared to 
experiment(see Figure 1 in Ref [13]). 

5. Conclusion 

We have preformed a unitary, full lattice QCD calculation of the r\ and tj' mass, and the isospin 
71 — 71 scattering length where the disconnected graph plays an important role. The r\ mass we 
obtain is in a good agreement with chiral pertubation theory. We still have large errors for the r\' 
mass, and the full I = % — % scattering length has such a big error that the result serves only 
as a bound on the magnitude . A clear signal can only be seen for small time separation in the 
disconnected diagram of % — % scattering. 

In the future, increasing the statistics alone could improve our results to some extent. Reducing 
the vacuum noise by doing the calculation on a larger volume or making the pion mass smaller 
could greatly improve the signal to noise ratio of the vacuum diagram. As more powerful computers 
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Figure 5: The K — K scattering length for 7 = 2 and 7 = 0. The solid line gives the tree level prediction. 

are under construction, we will make the calculation of the disconnected diagram more precise, and 
hope to attack the e' problem. 
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